Specific heat and entropy of TV-body nonextensive systems 



Hideo Hasegaws 

Department of Physics, Tokyo Gakugei University, Koganei, Tokyo 184-8501, Japan 

(Dated: August 13, 2010) 

Abstract 

We have studied finite iV-body L>-dimensional nonextensive ideal gases and harmonic oscillators, 
by using the maximum-entropy methods with the q- and normal averages (q: the entropic index). 
The validity range, specific heat and Tsallis entropy obtained by the two average methods are 
compared. Validity ranges of the q- and normal averages are < q < qu and q > qu, respectively, 
where qu = 1 + (nDN)^ 1 , q^ = 1 — (nDN + l)^ 1 and n = 1/2 (n = 1) for ideal gases (harmonic 
oscillators). The energy and specific heat in the q- and normal averages coincide with those in the 
Boltzmann-Gibbs statistics, although this coincidence does not hold for the fluctuation of energy. 
The Tsallis entropy for N\q — 1| S> 1 obtained by the g-average is quite different from that derived 
by the normal average, despite a fairly good agreement of the two results for \q — 1| <C 1- It has 
been pointed out that first-principles approaches previously proposed in the superstatistics yield 
additive iV-body entropy = NS^) which is in contrast with the nonadditive Tsallis entropy. 
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I. INTRODUCTION 



Many studies have been reported on the nonextensive thermodynamics since Tsallis pro- 
posed the so-call Tsallis entropy defined by [jJ-S] 



s q = J^- (J ( Pq y d n - 1) 



(i) 



where ks denotes the Boltzmann constant, p q the probability distribution function (PDF) 
and Q the phase space element. The Tsallis entropy is a one-parameter generalization of that 
in the Boltzmann-Gibbs (BG) statistics, to which Eq. ([!]) reduces in the limit of q = 1.0, 



Si = S B g = -k B / pi lnpi dVL. 



Entropy additivity is defined t 
independent systems A and B 



(2) 



rat an entropy S is additive if for any two probabilistically 



n 



p q (x 1 ,x 2 ) = Vq{xi)Vq{x2) for xi E A, x 2 E B, 



(3) 



we have 



S(A + B) = S(A) + {B). 



(4) 



7rom the definition, the BG entropy is additive whereas the Tsallis entropy is nonadditive 



4j because Eqs. (CQ) and (131) lead to 



S q (A + B) = S q (A) + S q (B) + (l-q)S q (A)S g (B). 



(5) 



For probabilistically independent, identical iV-body systems, an entropy S(N) is additive if 
we have S(N) = NS(1) for ViV. 

□ 

Entropy extensivity is a concept more subtle than the additivity in some sense [7|. An 
entropy S(N) of a given iV-body system is said to be extensive if 



S(N) 

< hm , T < oo, 



(6) 



i.e., if S(N) oc for ^> 1 [7]. Additivity depends only on the specific mathematical 
connection between the macroscopic entropy functional and PDFs of the system. Extensivity 
depends on this but also on the nature of the interaction and/or correlation of its elements. 
The distinction between additivity and extensivity has already been illustrated in simple 



probabilistic systems (see Fig. 1 in Ref. Q). For independent (or short-ranged interaction) 
systems, the BG entropy is extensive while the Tsallis entropy is nonextensive. In contrast, 
the BG entropy is nonextensive for systems with long-ranged interactions. The Tsallis 
entropy may be extensive for correlated systems such as the probabilistic system js], some 
fermionic system [9| , and Gaussian PDFs with correlation [10| . The nonextensivity is realized 
in four different classes of systems 4): (a) systems with long-range interactions, (b) small- 
scale systems with fluctuations of temperatures or energy dissipations, (c) systems coupled 
to finite-size thermal reservoir, and (d) multi-fractal systems. The Tsallis nonextensive 
statistics has been successfully applied to a wide range of nonextensive systems including 
physics, chemistry, astronomy, geophysics, biology, economics, and others 11]. 



The specific heat and Tsallis entropy in nonextensive systems such as ideal gases 



12|- 



181 ] and harmonic oscillators 18H20) have been calculated mostly by using the g-average 
which is employed in the un-normalized maximum-entropy method (MEM) [2], normalized 
MEM jsj, and the optimum Lagrange multiplier (OLM) MEM [6j. This is because the 
MEM with the g-average has been believed to be more superior than the original MEM 

1] with the normal average (4, 21 23| . Recently, however, it has been pointed out that 
thermodynamical quantities obtained by the g-average are unstable for a small change o 
the PDF, whereas those obtained by the normal average are stable 



26|]. Ref. [27] 



24 



has shown that expectation values may be self-consistently formulated with the form of 
normal averages in the generalized variational perturbation approximation. In contrast, it 
has been claimed that the Tsallis entropy and thermodynamical averages are robust for 
the g-average j^sj. The stability (robustness) of thermodynamical averages which is one of 
fundamental issues in the nonextensive thermodynamics, is currently controversial [29) . It 
is worthwhile to make calculations of thermodynamical quantities in iV-body nonextensive 
systems with the use of both the q- and normal averages. Although the N dependences of 
specific heats in ideal gases and harmonic oscillators have been calculated 12|- 20|], those of 
the Tsallis entropy have not been investigated, as far as we are aware of. The purpose of the 
present paper is twofold: (1) to calculate the N dependence of thermodynamical quantities 
of nonextensive ideal gases and harmonic oscillators and (2) to make a comparison between 
results with the q- and normal averages. Such calculations are ex pec ted to provide some 



insight to the current controversy on the two average methods 



24j-|29]. 



One of alternative approaches to nonextensive systems is the superstatistics 30|- [33] . 



3 



Complex nonextensive systems are expected to temporary and spatially fluctuate. In the su- 
perstatistics, it is assumed that locally the equilibrium state is described by the Boltzmann- 
Gibbs statistics, and that their global properties may be expressed by a superposition of 



them over some intensive physical quantity, e.g. the inverse temperature 



30|-|33|. Many 



applications of the concept of the superstatistics have been pointed out (for a recent re- 



view, see 



331] ). It is, however, not clear how to obtain the mixing probability distribution 



of fluctuating parameter from first principles. This problem is currently controversial and 



some attempts to this direction have been proposed [3J]-[37|]. Previous studies in the su- 
perstatistics have paid attention mostly to the stationary PDF, which generally has the 
non-Gaussian form. However, thermodynamical quantities such as energy and entropy have 
not been calculated in the superstatistics, which will be studied in this paper. 

The paper is organized as follows: In Sec. II, we briefly review the OLM-MEM with the 
g-average [6( and the original MEM with the normal average MEMs with the two 

average methods are applied to ideal gases and harmonic oscillators in Sees. Ill and IV, 
resp ectively. In Sec. V, the fluctuation of the energy and applications of the superstatistics 
[30[- 33j to iV-body systems are studied. Sec. VI is devoted to our conclusion. 



II. MAXIMUM-ENTROPY METHODS 
A. OLM-MEM with the g-average 

We consider an independent iV-body system whose hamiltonian is given by 

N 

H = (7) 

i=i 

In order to obtain the iV-variate density matrix p q from the Tsallis entropy with the OLM- 
MEM [6|, we impose the constraints given by 

1 = TV p„ (8) 

Cq 

with 

c g = Tr (p q y, (10) 
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where the bracket [-] q denotes the q average, Tr the trace and U q the g-averaged energy. The 



OLM-MEM yields the density matrix given by [6] 



pq = T q)l {H ~ Uq)] 



1/(1-?) 
+ 



with 



X q = Tr [1 - (1 - g) 7 (H - U q )\ 



1/(1-?) 



(11) 



(12) 



where [x] + = max(x, 0) and 7 expresses a Lagrange multiplier. Although /3 is conventionally 
used for the relevant Lagrange multiplier, we employ 7 to avoid a confusion with the inverse 
of physical temperature (T = l/fce/3). The Tsallis entropy defined by Eq. (jTJ is expressed 
in terms of c q by 



S„ 



q 



(13) 



where c q is given by Eq. f flQT) . 

Next we relate the Lagrange multiplier 7 to the physical temperature T. It is an open 
problem how to define the physical temperature in conformity with the zeroth law of ther- 
modynamics in the nonextensive statistics 



defining the physical temperature 



33|- 



43|. Among several schemes proposed for 



431 ]. the method of Ref. 



38 



this study. We assume that the physical temperature T is given by 
Appendix A, related discussions being also given in Sec. VB) 

1 1 OS, 



has been employed in 
(for details, see 



T 



C q dU q 



which yields 



7 



We may rewrite Eq. ( TlTl) as 



with 



Pq 



[l-(l-q)aH] 



1/(1-3) 



(14) 



(15) 



(16) 



Tr [1 



q)aH] 



1/(1-9) 



(17) 



U q is expressed by 



U g = —-Tr{[l-(l-q)aH)f 1 -' l) H}, 

with 

Vq = [l + (l-q)PU q ], 

where we adopt the relation: c q = (X q ) 1 ~ g = h / g (Z q ) 1 ^ q . 



(19) 



(20) 



B. Original MEM with the normal average 



n 



In the original MEM with the normal average [l|, we impose the constraints given by 
(tildes are attached to quantities relevant to the normal average) 



1 = Trpg, 
U q = (H) q = Tr {p q H}, 



(21) 
(22) 



where the bracket (■)„ expresses the normal average and U q the normal- averaged energy. 



The original MEM [1 



24| yields the density matrix give by 



l-{ q -l)^[H-U, 



1/(9-1) 



(23) 



with 



X q = Tr 



l-(q-ltf(H-U, 



1/(9-1) 



where 7 denotes a Lagrange multiplier. The Tsallis entropy is expressed by 



S a 



Cq ~ 1 

l-q- 



(24) 



(25) 



with 



c q = Tr (p g y. 



(26) 



We assume that the physical temperature T is given by [38| (for a detailed discussion 
see Appendix A) 

1 1 dS, 



T CqdU q 



qk B 7, 



(27) 
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which leads to 



7 







Equation (123]) may be rewritten as 



with 



1 



Z q = Tr 



q-1 



aH 



1/(9-1) 



q-1 



aH 



1/(9-1) 



(28) 



(29) 



(30) 



U q is expressed by 



a. 



l-{l/q-l)/3U q 



U q = -^Tr { 



I - I q - ) n // 

1 



1 V(9-l) 



H}. 



(31) 



(32) 



Calculation methods discussed in this section will be applied to classical ideal gases and 
harmonic oscillators in Sees. Ill and IV, respectively. 

III. IDEAL GASES 
A. g-average 



We consider hamiltonian for iV-body .D- dimensional ideal gases given by 

N N 



i=l 



i=l 



2m' 



(33) 



where m stands for the mass and p { momentum of ideal gases. By using the exact method 



13 



481 ] . we obtain the partition function Z q in Eq. ( ITT]) (for details, see Appendix Bl), 



yN 


2-7T m 


DN/2 r . 


1 DN \ 
q-1 2 ) 


M h DN 


_(q-l)a 






yN 






N\ h DN 

yN 


\ a ) 
2ivm 


DN/2 




N\h DN 


(l-q)a 


r( 





for 1 < q < 3, 

for 5 = 1, 
for q < 1, 



(34) 
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where h stands for the Planck constant and V a volume of the system. The average energy 
in Eq. ffl9l) is given by (see Appendix Bl) 

DN 



[DN(l-q) + 2]a 



From Eqs. ffl8|) and f )35|) . U q and a are self-consistently determined as 



U„ 



a 



DN _ k B DNT 
~W ~ 2 ' 



1 + 



DN , 



l-q) 



A conceivable q value is given by 



< q < 1 + 



(35) 

(36) 
(37) 

(38) 



because p q given by Eq. f fl6|) has the probability properties with a > and because the 
Tsallis entropy is stable for q > 2l|. Equation (|36|) leads to the specific heat given by 

k B DN 



a, 



(39) 



which agrees with that in the BG statistics. The specific heat of nonextensive ideal gases 
was discussed in Ref. [l^] by using the g-average with the un- normalized MEM where 
a possibility of the negative specific heat was pointed out. Our result of Eq. (139]) is in 
agreement with previous studies using the OLM-MEM Q-QQ 19] . 
The Tsallis entropy may be calculated by 



S q = k B ' 



with 



Q 



DN , 



The BG entropy in the limit of q = 1.0 is given by 

S 1 = NS[ 1] , 



with 



(i) 



27rm 



In 



a 



D 

i + 1 



(a = V/N). 



(40) 

(41) 
(42) 

(43) 

(44) 
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By using the formula: 



In 



T{z + a) 
T(z)z a 



a(l — a) 

2z 



for \z\ — >• oo, 



(45) 



we obtain S q for \q — 1| -C 1.0 given by 



5, = NSl 1] - (q - 1 



iV 2 DiV 



+ 



(46) 



where the (g — 1) term includes the 0(N 2 ) contribution showing the nonadditivity. 
The Tsallis entropy for large N(l — q) is given by 



q ~ 



for N(l - q) > 1, 



(47) 



where we adopt the relation: lnT^) ~ (z — 1/2) haz — z + (1/2) ln27r + ■■ for \z\ — > 00. 



B. Normal average 

For ideal gases whose hamiltonian is given by Eq. (|33|) . we obtain (for details, see 
Appendix B2) 



yN 


2n m 


DN/2 r( _i_ +l) 


N\ h»*> 


_(q-l)a_ 


r(^r+^+i) 






N\ h DN 

yN 


\ a ) 
2irm 




N\ h DN 

V 


_(l-q)a_ 





for 1 < q < 3, 
for q — 1, 
for q < 1. 



U q in Eq. (|32|) may be expressed in terms of a (see Appendix B2) 

u. 



[DN(q-l) + 2q]a 
From Eqs. (13T1) and (j4"9l . £7q and a are self-consistently determined as 

DN k B DNT 



a 



2/3 2 

ff+¥(g-i)' 



A conceivable q value is given by 

q > 1 



(¥ + 1) 



<?L = — , 



(48) 



(49) 



(50) 
(51) 



(52) 
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where qu is given by Eq. (|38|) . Equation ( |50|) yields specific heat given by 

^ k B DN 

which is the same as that of the BG statistics. 
The Tsallis entropy is expressed by 



S n — k 



B 



with 



Cg-1 

l-q 



{Zq) 1 -* 



v„ 



DN 



Vn 



!--)+! = vi/ q , 



where v q is given by Eq. fl42|) . 

For \q — 1| <C 1.0, we obtain S q given by 

S, = NS™ - (q - 1) 



+ -, 



(53) 



(54) 



(55) 
(56) 



(57) 



where is given by Eq. f )4*4|) . Equation ( 15 7p is the same as Eq. ( H6|) except for the O(DN) 
term in the bracket. 

For large iV(g — 1), S q is given by 



Sq — 



k 



B 



q-l 



2q e -{ q -i)NS<p /k B 
N(q-l) 



1 



for N{q - 1) > 1, 



(58) 



which is quite different from S g in Eq. fll7|) obtained by the g-average. 



C. Model calculations 

We will present some numerical calculations of the Tsallis entropy of ideal gases with 
D — 1. Equations ( 138]) and (1521) show that the conceivable q values for q- and normal averages 
are < q < qu and q > qu, respectively, where qu = 1 + 2/DN and qu = l — 2/(DN + 2) (see 
the inset of Fig. [T]). qu and qu approach unity for N — > oo. Although the iV dependence 
of the Tsallis entropy in the g-average is rather different from that in the normal average, 
as Eqs. (14T|) and (158]) show, both results are in fairly good agreement for q L < q < qu 
where both the q- and normal averages are valid. Dashed and solid curves in Fig. [1] show 

10 



the entropies calculated by the g- and normal averages, respectively, for ideal gases with 
N = 100 and T/T = 1.0 where T = h 2 /k B m (q L = 0.9804 and q v = 1.020). The difference 
between the two entropies is small for q L < q < q v as shown by the bold curve in Fig. [T] 
except for q < q L where S g is divergently increased. 

Figure [2] shows the three-dimensional plot of Sq/ksN calculated by the g-average as 
functions of q and N for T/Tq = 1.0. With increasing N, S q is considerably increased, as 
Eq. (j47p shows: note the logarithmic scale in the ordinate. 

The temperature dependence of S q calculated by the g-average is shown in Figs. [3] and 
HJ in the former Sg/kgN for N = 100 is plotted as functions of T and q, while in the latter 
Sg/ksN for q = 0.8 is plotted as functions of T and N. With increasing the temperature, 
the Tsallis entropy is increased, as expected. This temperature dependence becomes more 
significant for larger N. 

Figure [5] shows S q /k B N calculated by the normal average for T/T = 1.0 as functions 
of q and N. The magnitude of Sq/ksN is significantly decreased with increasing N and/or 
q, as Eq. (1581) shows. This behavior is quite different from that obtained by the g-average 
shown in Fig. [3 

Figure [6] shows S q / kgN calculated by the normal average as functions of T and q, whereas 
Fig. [7| shows Sq/ksN as functions of T and N. With increasing the temperature, S q is 
increased. The temperature dependence of S q becomes weaker for larger N and/or larger q. 



IV. HARMONIC OSCILLATORS 



1. q-average 



We consider hamiltonian for iV-body D- dimensional harmonic oscillators given by 

N N / 2 2 2\ 

i=l i=l v 7 

where m, u, r» and p { express the mass, oscillator frequency, position and momentum, 
respectively. From Eq. (JIT]) , we obtain (see Appendix Bl) 



h DN 



2- 



(q-l)uia 



_J_ (2e.\ 

h DN \wa) 



DN 



DN r (_i__ DJV ) 
n^t) 



1 



2tt 



(1— q)uia 



DN 



. 1-9 



r (T L - q + 1+DN ) 
11 



for 1 < q < 3, 

for 5 = 1, 
for g < 1. 



(60) 



U q in Eq. (fT9l) is expressed in terms of a (see Appendix Bl), 

DN 



U„ 



q ' [DN(l-q) + l]a' 
From Eqs. (fT8|) and ( IBTj) . {7 g and a are self-consistently determined as 

DN 

U n = —r- = k B DNT, 







a 



ft 



1 + DN(1 - q) 
Equation (|62|) yields the specific heat given by 



for < q < 1 + = q v . 



C q = k B DN. 



(61) 

(62) 
(63) 

(64) 



The specific heat of harmonic oscillators was discussed in Ref. 19( with the use of un- 
normalized J5] and normalized MEMs |2j. Their expressions for the internal energy U q have 
rather complicated q dependence [see Eqs. (12) and (22) of Ref. [19|]. Equation (151")) 
expresses the Dulong-Petit specific heat: Cbg = k B DN in the BG statistics, which is in 
agreement with Ref. 18]. 

The Tsallis entropy is given by Eq. (j4"0]) with 



Uq = DN(1 -q) + l. 



The BG entropy in the limit of q = 1.0 is given by 



Si = NS\ 



(i) 



with 



(i) 



k B D 



In 



(t^) +1 



(65) 
(66) 

(67) 

(68) 



2. Normal average 



For harmonic oscillators under consideration, we obtain (see Appendix B2) 



h DN 



2- 



(q-l)uia 



DN 



,9-1 



J_ ( 2JL ) 



DN 



h DN 



2tt 



(1— q)uia 



V(^+DN+1) 



DN r (_l__ DA r) 



for 1 < q < 3, 

for 5 = 1, 
for q < 1. 



(69) 
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We may express U q given by Eq. ( 132|) in terms of a as given by (see Appendix B2) 

DN 

Uq = [DN(q-l) + q]a (70) 
From Eqs. (l3Tj) and (1701) . f/ g and a are self-consistently determined as 

DN 

U q = — = k B DNT, (71) 

" = q + DN(q-l) for <? > 1 " P+i) = 9l = 1/qu, (72) 
where ^ is given by Eq. ( 1631) . We obtain the specific heat given by 

C q = k B DN. (73) 

Calculated specific heat is the same as the Boltzmann-Gibbs result of Cbg — k B DN. 
The Tsallis entropy may be calculated by Eq. (|54p with 



c, 



V 1 



(74) 



v q = DN (l - + 1 = u 1/g , (75) 

where v q is given by Eq. (|66|) . 

The properties of the specific heat and Tsallis entropies of harmonic oscillators are es- 
sentially the same as those of ideal gases if we read DN/2 — >■ DN and m — > 1/u. 

V. DISCUSSION 

A. Fluctuation of the energy 

It is interesting to examine the fluctuation of energy in nonextensive systems, which has 



been studied in Refs. 



18 



44 



451 ]. The g-average of H 2 is given by (see Appendix Bl) 



.., 7]DN(r]DN + 1)[1 + T]DN(l-q)\ 



[H]q F[l + ( V DN+l)(l-q)} ' (76) 

where rj — 1/2 and rj = 1 for ideal gases and harmonic oscillators, respectively. The relative 
fluctuation of the energy is given by 



AU, _ JVP-V$, 

U q U q ^/r)DN[l + (rjDN + 1)(1 - q)] ' 

f i for q = 1.0, 
^ for DN(1 - q) » 1. 
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Equation (1771) agrees with the previous results with the OLM-MEM Q, fl |45|. Equation 



f8j) shows that when the nonextensivity exists, the relative fluctuation is proportional to 
1/N for large N instead of 1/y/N in the BG statistics. 
The normal average of H 2 is given by (see Appendix B2) 

V DN (rjDN + !)[! + (yDN + l)(g - 1)] 
^ ;<? /3 2 [l + (77L>iV + 2)(g-l)] ' 1 ' 

which leads to the relative fluctuation of the energy given by 



AU q \/{H 2 -U*) q x 



U q U q VvDN [l + ( V DN + l)(l-l/q)] 



(80) 



oc fOT ^ = L °' (81) 

forDJV(g-l)»l. 

Equation ( 180]) has the similar DN dependence to Eq. (ITT)) . We note that Eqs. ( ITTj) and 
( IHUj) have the reciprocal symmetry: g -H- 1/g. 

Figure [8] shows the 3d plot of AU q /U q as functions of q and (rjDN) calculated by the q 



average. The DN dependence of AU q /U q for large N is changed from 1/yDN to 1/DN 
when the nonextensivity is introduced as Eq. (ITBl) shows. We notice also that the magnitude 
of AU q /U q is much decreased by the nonextensivity. The DN dependence of the relative 
fluctuation of energy in the normal average is given in Fig. |8] if we read q — > 1/q. 



B. Superstatistics 



1. Probability distribution function 



In Sec. II, we have derived the PDF by using the MEMs. It is instructive to obtain the 



PDF in the superstatistics 



30]- 



331 ] . The PDF of .D-dimensional iV-body ideal gases in the 



locally equilibrium state with a given inverse temperature (3 is expressed by 

AT! /h^\ DN ' 2 



n m (p;« = ^( 



2irm J 



-/9E£iP?/2m 



(82) 



where p = {p { } (i = 1 to N). We assume that the distribution of (3 is given by the 
^-distribution with the rank n 30|, l3l| . 

/(/?) 



r(n/2) \2f3 



_^ | on/2-1 e -n/3/2/3^ 



n/2 



(83) 
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where (3 stands for the average of over f((3): (f3)p = (3. The PDF averaged over f((3) is 
given by 



POO 

p( N \p) = / nW(p;/3)/(/3)d/J, 
Jo 



1 

Z 



V n J ^-^ 2m 

7 i=i 



-(DN+n)/2 



(84) 
(85) 



with 



V N ( nmn 



DN/2 



r(n/2) 



(86) 



iV! V /i 2 /? 7 V{DN/2 + n/2) ' 
The PDF of harmonic oscillators in the local equilibrium state with a given inverse 
temperature (3 is expressed by 



hu/3 
~2^ 



DN 



e -^E£i(P?/2m+ma; 2 r^/2)_ 



The PDF averaged over /(/?) is given by 



P iN) (p,x) 



nW(p,x;P)f(P)dp, 



1 

z 



N 



2(3\ ^ Z' p 2 mu 2 r 2 



E 



" / z — ' V 2m ~^ 2 



-(DN+n/2) 



(87) 

(88) 
(89) 



with 



r 



~// \ DN r(n/2) 



(90) 



Vmv r(L>i\r + n/2) ' 

It is easy to see that PDFs given by Eqs. ( 185]) and ( 1891) have the same structure as those 
given by Eqs. (fTBj) and (1291 derived by the MEMs. For example, Eq. ( 1551 for ideal gases 
may be rewritten as 



p iN \p) oc 



1 - (1 - q) 



DN + n 



n 



\ N r> 2 

r ^2m 

' i=i 



1/(1-9) 



(91) 



with 



q = 1 + 



DN + n 



(92) 



or equivalent ly 



p (iV) (p) oc 



g- l\ [DN + n- 2 



n 



N 



^E 



i=i 



Pt_ 

2m 



i/(g-i) 



(93) 
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with 

9 = 1 - • ( 94 ) 

DN+n y ' 

Equations (|9Tj) and ( 193]) are similar to Eqs. (|T6|) and ( 129|) . respectively, although the ranges 

of conceivable g values are not necessarily the same. 



2. Energy and entropy 



The ^-distribution in Eq. (1831) is given ad hoc 30]-[33|. Several first-principles ap 



proaches have been proposed to determine the optimum distribution of (3 in the su- 
perstatistics j34]-j33|. Ref. considered the local energy E((3) = (H) n and entropy 
S((3) = — (lnn) n for local equilibrium states with an inverse temperature (3, and then 
obtained the energy and entropy of the system given by E = (E((3))p and S = {S((3))p, 
where (-)n and (-)p express averages over IT(/?) and f((3), respectively. Equations (I3B"|) and 
( |44l) for ideal gases lead to 



SiP) 



DN 

k B N 



In 



2irm 
~Wf3 



In cr + 



D 



+ 1 



from which we obtain 
DN 



E 

S 



2(3 
k B N 



n 



n 



D\ (2-nm 



+ In ex + ( — + 1 



2~ 



k B DN 



(95) 
(96) 

(97) 
, (98) 



ip{x) standing for the poli-gamma function. Similarly from Eqs. f )62|) and f )68|) for harmonic 
oscillators, we obtain 



SO?) 



DN 
k B DN 



ln 



2tt 



and then 



5 



In 



2vr 
fed? 



, n\ (n 



(99) 
(100) 

(101) 
(102) 
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We note that results given by Eqs. (1571) . (!98|) . f llOip and (I102p are different from their 
counterparts obtained by the MEMs. The similar result is obtained in the maximum-entropy 
approach to the superstatistics in Ref. 36] based on the conditional entropy. When we 
employ, in place of the gamma distribution of Eq. (153"j) . the inverse-gamma distribution 
proposed in Ref. [46j , 

n/2 



/t(/3) 







r(n/2) 



n/3 
2 



-n/2-2 -n/3/2/3 



(103) 



we obtain for ideal gases, 
L>iV /n + 2 s 



5 



2^ 



n 



27rm 



+ ln<r 



D 

2 +1 



, n 



ln[™ 



(104) 
(105) 



We note that entropies in Eqs. (l98lh fl 1 2 [) and fl 1 5 j) derived by the maximum-entropy ap- 

371 ] are additive, which is in contrast with the nonadditive 



36 



proaches to the superstatistics 
Tsallis entropy defined by Eq. (p]). This result is not modified even when we adopt any dis- 
tribution function of j3, because S(/3) is given by S(/3) = ksNh^fi) [see Eqs. ( 1961) and ( HOOD 
and then S = ksN '(h(/3)) p, where h((3) depends on /3 but on N. The superstatistics 3oT 31 
was introduced to account for the non-Gaussian distribution which is well described by the 
g-Gaussian (g-exponential) form originating from the nonadditive Tsallis entropy [l). We 
then naturally expect that the entropy obtained by the superstatistics is also nonadditive. 
It is, however, not the case in the first-principles approaches to superstatistics jsij Q as 
mentioned above. 

Alternatively, we may calculate the energy and entropy by using the partition function 
Z given by Eqs. (186"]) or (|90|) . For ideal gases, we obtain 

DN 



E' 
S' 



_dZ_ 

~~d(3 ~ ~2/T' 
k B QnZ + pE') 

k B N (?) h, 



(106) 



2irm 
1^3 



, D 

+ lna+ — 



+ k 



B 



DN 



n 



In - + In T - 



n 



lnT 



DN n 
~2~ + 2 



(107) 



E' agrees with Eq. ( 1361) or (|50|) if we read (3 = (3. Although S' is nonadditive, its 
dependence is different from S q derived by the MEM in Sec. III. 
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C. The physical temperature 



As mentioned in the introduction, there are several proposals in relating the Lagrange 



multiplier to the physical temperature even within the OLM-MEM 
we have adopted the scheme proposed in Ref. 



43j . In this paper, 



. One of the alternatives is to define the 



physical temperature by 1/T = dS q /dU q (3(, yielding for ideal gases (with the OLM-MEM), 



1 

f 



k B J c, 



q-> 



(10$ 



which is different from Eq. (fl4"l) by a factor of c q . Then the Lagrange multiplier is determined 

as 

1 



7 



and the energy is given by 



k B T c q 



DNk B T c q 



(109) 



(110) 



The expression for the specific heat becomes complex because c q in Eq. (IllOp is temperature 
dependent: Z q in c q has the temperature dependence through 7. 



VI. CONCLUDING REMARKS 



Specific heats and entropies in nonextensive ideal gases and harmonic oscillators have been 
calculated by using the MEMs with the q- and normal averages and by the superstatistics. 
Obtained results of ideal gases are summarized in Table 1: those of harmonic oscillators are 
similarly given if we read DN/2 — > DN. Our calculations have shown the followings: 

(i) The MEMs with the q- and normal averages are valid for < q < qu and qi < q < 00, 
respectively, where qu = 1 + 1/r/DN and q^ = \ — \j{r\DN + 1), 

(ii) Energy and specific heat obtained by the q average coincide with those derived by the 
normal average and in the BG statistics, 

(iii) the Tsallis entropy obtained by the q- and normal averages are in fairly good agreement 
for q L <q < qu, 

(iv) the Tsallis entropy in the g-average shows the exponential N dependence: S q /N ~ 
e ('i--q)NS 1 ) f Qr — ^> 1, which is quite different from that in the normal average: 
S q /N ~ l/(q - 1)N for N(q - 1) > 1, and 
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TABLE I: Specific heat C q and the Tsallis entropy S q of ideal gases [S^ is given 



by 



36 



/qs. 



m the x 2 



qu = 1 + 2/DN and g L = 1 - l/(DN/2 + 1) = l/gjy. In the superstatistics 
distribution with rank n is employed: G(n) = ip(n/2) — ln(n/2) oc 1/n for n S> 1. -/Wl: not 
applicable] 



method 


range 




S q (for \q - 1| < 1.0) 


S g (for A%-1| > 1) 


g-average 


< q < qu 


k B DN 
2 


AT^-^-l) 






k B N (l-q)NSP Iks 
2 e 


normal average 


q L < q < oo 


k B DN 
2 


ATSf )-(<,-!) 




2 + ^" 


(9-1) 


Superstatistics 


n > 2 


k B DNn 
2(n-2) 




JV \s{ 1] kB 2 D G(n)] 



(v) the entropy in first-principles approaches to the superstatistics of Refs. [36|, |37| is additive 
(gCAQ _ jysi 1 )^ while the Tsallis entropy is nonadditive. 

Items (i), (iii) and (iv) are consistent with results obtained with the g-Gaussian PDF \^\ . 
The item (ii) shows that the rule of energy equipartition holds in both BG and nonextensive 
statistics. The item (iii) implies that the calculated properties of the entropy considerably 
depend on the adopted average methods 0|. It would be necessary to further examine 



q- and normal averages from various viewpoints 26|]- 29] and to develop a first-principles 
superstatistical method which yields an entropy in conformity with the Tsallis one. 
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Appendix: A. Derivation of the physical temperature 



1. (/-average 



We will derive the physical temperature, defining 



X q = Tr [l-(l- g ) 7 (tf-[/ 9 )] 



V(l-9) 
+ 



Tr [l-(l- q ) 7 (H-U q 



19/(1 



where Tr denotes the trace. We obtain 



1-9 



The physical temperature T is assumed to be given by 

1 1 8S n 



T c q dU q 

k B (dX t 



X q \dU q 



yielding 



7 



where we employ the relation: dX q /dU q = 7 X q . 



2. Normal average 



For the normal average, we define X q and c„ by 



X q = Tr 



l-(q-ltf(H-U< 



Tr 



1/(9-1) 



9/(- 



A simple calculation leads to 
1 



1 - (9 — 1)7 [H —U, 



^Tr 

^Tr{[l-(g-l) 7 (H-U Q 



9/(9-1) 



1/(9-1) 



l-(q 
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where we employ the definition of U q given by 

U q = ^Tr{[l-(q-l)^H]f^H}. 

Equation (1A8|) leads to 

c q = {x q y-<i = {v q )-\z q Y-«. 

With the use of Eqs. (IA10I) and (Elj), the physical temperature is given 

1 1 dS q 



t c Q du: 



q UU q 

k B \ dX n 



x q l dU q 



leading to 



~ 
7 = - 
Q 



where we employ the relation: dX q /dU q = q 7 X q . 

Appendix: B. Evaluations of averages by the exact approach 
1. The (/-average 

We first discuss evaluations by the g-average given by 

Z q (a) = Tr [l-{l-q)aH} l ^- q \ 

Q q (a) = —L-T T {{I - (I - q) a Hf^ Q}, 
v q A q [a) 



by using the exact expressions for the gamma function 

1 



13 



48|: 



CO 



y- s = —— I u s - l e~ yu du for s > 0, 

i(s) 







i 
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T(s + 1) f {-t)- s ~ l e~ yt dt for s > 0, 
Jc 
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where Q denotes an arbitrary operator and C the Hankel path in the complex plane. We 



obtain I13L 



48| 



Z g (a) 



f^u^^-^e^Z^q - l)au] du for q > 1.0, 

+ l) / (-t)-V(i-9)-i e -%[_(l - q)at] dt for q < 1.0, 



(B5) 



Qq(0l) 



7 , ! w 5 s f^u^-^e-vZtiiq - l)au] Q x [{q - l)au] du for q > 1.0, 



2-KVq Z q (a) 



r + l) J c (-i)-^ 1 -^- 1 e- t Z l [-(l - q)at] Q![-(l - q)at] dt, 



for q < 1.0, 
(B6) 



where 



Z\{a) = Tie 

QM) 



-aH 



Tv{e- aH Q} 



(B7) 
(B8) 



and v q is given by Eq. (|42|) and (|66|) for ideal gases and harmonic oscillators, respectively. 



2. The normal average 



Next we discuss evaluations by the normal average given by 



Z q (a) = Tr [1 - (q- 1)5 H] 



1/(9-1) 



— Tr{[l - (q- l)aH} 1/{q - 1] Q} 



By using Eqs. (1B3|) and flB4j) . we obtain 



(B9) 
(BIO) 



+ l) / (_t)-V(9-i)-i e -t - 1)5*] rft for g > 1.0, 

j-« u i/(i-,)-i e -« _ g)£m] du for ? < i.o, 



(Bll) 
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»r[i/( g -i)+i] 



gtli J c (_ t )-i/{q-i)-i e -t Z x \-{q - l)at] Qi[-(q - l)dt] dt for g > 1.0, 



z q ^n m - q)] /"ttva-ffj-ig-i. zr ; [U - g)a«] gx[(i - g )st] rf M for 9 < 1.0, 

(B12) 



where Zx(a) and ?7i(a) are given by Eqs. (1B7P and (1B8j) . respectively. 
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FIG. 1: (Color online) Sq/ksN of ideal gases as a function of q for N = 100 calculated by the 
q average (g-av.: the dashed curve) and the normal average (N-av.: the solid curve), and the 
difference between the two entropies (the bold curve) (D = 1, T/Tq = 1.0: To = h 2 /fern). The 
inset shows the validity range of the two averages in the q-N space: the MEMs with the q- and 
normal averages are valid for < q < qjj and q > qL, respectively, and for qL < Q < Qu both the 
averages are valid: the dashed line denotes N = 100. 

FIG. 2: (Color online) Sq/ksN as functions of q and N obtained by the q- average (-0 = 1, 
T/Tq = 1.0): note the logarithmic scale in the ordinate and abscissa for N. 

FIG. 4: (Color online) Sq/ksN as functions of T and N obtained by the g-average (D = 1, 
q = 0.8): note logarithmic scales in the ordinate and abscissa for N. 



FIG. 5: (Color online) Sq/ksN as functions of q- and N obtained by the normal average (D = 1, 
T/Tq = 1.0): note the logarithmic scale in the abscissa for N. 



FIG. 6: (Color online) Sg/ksN as functions of q and T obtained by the normal average (D = 1, 
N = 100). 



FIG. 3: (Color online) Sq/ksN as functions of q and T obtained by the g-average (D = 1, 
N = 100): note the logarithmic scale in the ordinate. 
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FIG. 7: (Color online) S q /ksN as functions of T and N obtained by the normal average (D = 1, 
q = 1.2): note the logarithmic scale in the abscissa for N. 



FIG. 8: (Color online) The relative fluctuation of AU q /U q as functions of q and (r]DN) calculated 
by the q average: rj = 1/2 (r/ = 1) for ideal gases (harmonic oscillators). 
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